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Abstract 

We study the number of connected spanning subgraphs fd,b{ n ) on the generalized Sierpinski 
gasket SGd,b(n) at stage n with dimension d equal to two, three and four for 6 = 2, and layer 
' b equal to three and four for d = 2. The upper and lower bounds for the asymptotic growth 

constant, defined as zsc db = hm^-nx) ^fd,b{n)/v where v is the number of vertices, on SG^^n) 
with 6 = 2, 3, 4 are derived in terms of the results at a certain stage. The numerical values of zsG d b 
are obtained. 
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I. INTRODUCTION 



The enumeration of the number of connected spanning subgraphs Ncssg(G) on a graph 
G is a problem of interest in mathematics It is well known that the number of con- 
nected subgraphs is given by the Tutte polynomial T(G,x,y) evaluated at x = 1, y = 2 
{2I. Alternatively, it corresponds to the partition function of the g-state Potts model in 
statistical mechanics with q = and the temperature variable v = e^ J — 1 = 1, where J 
is the spin-spin coupling and j3 = (ksT) -1 . Some recent studies on the enumeration of 



connected spanning subgraphs and the calcu 



ation of their asymptotic growth constants on 



regular lattices were carried out in Refs. 0, \a , 0, [f], 0] . It is of interest to consider connected 
spanning subgraphs on self-similar fractal lattices which have scaling invariance rather than 
translational invariance. Fractals are geometric structures of (generally noninteger) Haus- 
dorff dimension realized by repeated construction of an elementary shape on progressively 
smaller length scales . A well-known example of a fractal is the Sierpinski gasket. We 
shall derive the recursion relations for the numbers of connected spanning subgraphs on the 
Sierpinski gasket with dimension equal to two, three and four, and determine the asymptotic 
growth constants. We shall also consider the number of connected spanning subgraphs on 
the generalized Sierpinski gasket with dimension equal to two. 

II. PRELIMINARIES 

We first recall some relevant definitions for connected spanning subgraphs and the Sier- 
pinski gasket in this section. A connected graph (without loops) G = (V, E) is defined by 



its vertex (site) and edge (bond) sets V and E |l0|, LL3J. Let v(G) = \V\ be the number of 
vertices and e(G) — \E\ the number of edges in G. A spanning subgraph G' is a subgraph 
of G with the same vertex set V and an edge set E' C E. A connected spanning subgraph 
on G is a spanning subgraph of G that remains connected. In general, there can be cycles 
in a connected spanning subgraph. It is called a spanning tree when there is no cycles. The 
degree or coordination number fcj of a vertex i>j e V is the number of edges attached to it. 
A A;-regular graph is a graph with the property that each of its vertices has the same degree 
k. In general, one can associate an edge weight Xij to each edge connecting adjacent vertices 
Vi and Vj. For simplicity, all edge weights are set to one throughout this paper. 
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When the number of connected spanning subgraphs Ncssg(G) grows exponentially with 
v{G) as v{G) — > oo, there exists a constant zq describing this exponential growth: 

\nN CSSG (G) 



zg 



lim 

r(G)-+oo 



v(G) 



(2.1) 



where G, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit. 

The construction of the two-dimensional Sierpinski gasket SG2{n) at stage n is shown 
in Fig. CD At stage n — 0, it is an equilateral triangle; while stage (n + 1) is obtained by 
the juxtaposition of three n-stage structures. In general, the Sierpinski gaskets SGd can be 
built in any Euclidean dimension d with fractal dimension D = ln(d + l)/ln2 12J. For the 
Sierpinski gasket SGd(n), the numbers of edges and vertices are given by 



e(SG d (n)) 



d 



(d + l)* = id + l) 



n+l 



v{SG d {n)) 



d+1 



[(d+ir+i] 



(2.2) 



(2.3) 



Except the + outmost vertices which have degree d, all other vertices of SGd(n) have 
degree 2d. In the large n limit, SGd is 2d-regular. 




A 

5G 2 (0) 5G 2 (1) 5G 2 (2) 




5G 2 (3) 



FIG. 1: The first four stages n = 0, 1,2, 3 of the two-dimensional Sierpinski gasket S'G^n). 



The Sierpinski gasket can be generalized, denoted by SG d ,b(n), by introducing the side 
length b which is an integer larger or equal to two [13j . The generalized Sierpinski gasket at 
stage (n + l) is constructed from b layers of stage n hypertetrahedrons. The two-dimensional 
S , G 2j fe(n) with b = 3 at stage n = 1, 2 and b = 4 at stage n = 1 are illustrated in Fig. [21 
The ordinary Sierpinski gasket SGd{n) corresponds to the 6 = 2 case, where the index b is 
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5G 2)3 (1) SG 2i3 (2) SG 2A (l) 

FIG. 2: The generalized two-dimensional Sicrpinski gasket SG2,b{n) with b = 3 at stage n = 1, 2 and 6 = 4 
at stage n = 1. 

neglected for simplicity. The Hausdorff dimension for SGd.b is given by D = In ( b+d d ~ 1 ) /In 6 



13j| . Notice that SGd,b is not fc-regular even in the thermodynamic limit. 



III. THE NUMBER OF CONNECTED SPANNING SUBGRAPHS ON SG 2 {n) 

In this section we derive the asymptotic growth constant for the number of connected 
spanning subgraphs on the two-dimensional Sierpinski gasket SG2{n) in detail. Let us start 
with the definitions of the quantities to be used. 

Definition III.l Consider the generalized two-dimensional Sierpinski gasket SG2,b{n) at 
stage n. (i) Define f2,b{ n ) = NcssG{SG2 t b{n)) o,s the number of connected spanning sub- 
graphs, (ii) Define g2,b(n) as the number of spanning subgraphs with two connected com- 
ponents such that one certain outmost vertex, say the topmost vertex as illustrated in Fig. 
[3] /or ordinary Sierpinski gasket, belongs to one component and the other two outmost ver- 
tices belong to another component. (Hi) Define h2,b(n) as the number of spanning subgraphs 
with three connected components such that each of the outmost vertices belongs to a different 
component. 

Since we only consider the ordinary Sierpinski gasket in this section, we use the notations 
/2 ( n ) > 92 (n) and h 2 (n) for simplicity. They are illustrated in Fig. [31 where only the outmost 
vertices are shown. Because of rotational symmetry, there are three possible g2{n). The 
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initial values at stage zero are /2(C)) = 4, #2(0) = 1 and /12(C)) = 1. The purpose of this 
section is to obtain the asymptotic behavior of f 2 (n) as follows. The three quantities f^fo), 
g 2 {n) and h 2 (n) satisfy recursion relations. 




/2O) 9z{n) h 2 {n) 

FIG. 3: Illustration for the connected spanning subgraphs f2(n), 92{'n) and h%(rt). The two outmost vertices 
at the ends of a solid line belong to one component, while the two outmost vertices at the ends of a dot line 
belong to separated components. 



Lemma III.l For any non-negative integer n, 

f 2 (n + l) = f*(n) + 6f>(n)g 2 (n) , (3.1) 

g 2 (n + 1) = fl{n)g 2 {n) + f 2 {n)h 2 {n) + 7f 2 (n)g 2 2 (n) , (3.2) 

h 2 (n + 1) = 3f 2 (n)g 2 2 (n) + 12 f 2 (n)g 2 (n)h 2 (n) + Ug 3 2 (n) . (3.3) 

Proof The Sierpinski gasket SG 2 (n + 1) is composed of three SG 2 (n) with three pairs 
of vertices identified. The number f 2 (n + 1) consists of one configuration where all three 
SG 2 (n) belong to the class that is enumerated by f 2 (n), and six configurations where one 
of the SG 2 (n) belongs to the class enumerated by g 2 (n) and the other two belong to the 
class enumerated by f 2 {n) as illustrated in Fig. HI Eq. (13. ip is verified by adding these 
configurations. 

FIG. 4: Illustration for the expression of f^n + 1). The multiplication of three on the right-hand-side 
corresponds to the three possible orientations of SG2(n +1). 
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LI = AA + AA + M + AS 2+ AA x2+ M 2 

FIG. 5: Illustration for the expression of 52 (»i + !)■ The multiplication of two on the right-hand-side 
corresponds to the reflection symmetry with respect to the central vertical axis. 




FIG. 6: Illustration for the expression of fi2(n + 1). The multiplication of three on the right-hand-side 
corresponds to the three possible orientations SG2(n + 1). 

Similarly, g 2 {n + 1) and h 2 (n + 1) for SG2(n + 1) can be obtained with appropriate 
configurations of its three constituting SG2{n) as illustrated in Figs. OandOto verify Eqs. 
(13. 2p and (13. 3p . respectively. 



□ 

The values of /2(H), g2(n), h2(n) for small n can be evaluated recursively by Eqs. ( 13. ip . 
A3. 21) . (13.31) as listed in Table [B These numbers grow exponentially, and do not have simple 
integer factorizations, in contrast to the corresponding results for the number of spanning 
trees [ijj]. To estimate the value of the asymptotic growth constant defined in Eq. (12. ip . we 
need the following lemmas. For the generalized two-dimensional Sierpinski gasket SG2,b{n), 
define the ratios 

Ol2,b{n) = p-r , P2,b{n) = - — p-r , (3.4) 

92,b{n) h 2 ,b{n) 

where n is a non-negative integer. For the ordinary Sierpinski gasket in this section, they 
are simplified to be ct2 (n) , fa {n) and their values for small n are listed in Table HT1 



Lemma III. 2 For any n > 0, 

3/3 2 (n) < a 2 (n) < A(3 2 (n) . (3.5) 
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TABLE I: The first few values of /^(n), <7 2 (n), h2(n). 



n 





1 


2 


3 


/ 2 (n) 


4 


160 


13,312,000 


10,293,452,839,321,600,000,000 


92 (n) 


1 


60 


7,462,400 


8,864,355,990,896,640,000,000 


ha(n) 


1 


74 


13,276,800 


23,868,720,258,482,176,000,000 



TABLE II: The first few values of a 2 (n), P%{n)- The last digits given are rounded off. 



n 





1 


2 


3 


4 


a 2 (n) 


4 


2.66666666666667 


1.78387650085763 


1.16121835019855 


0.736689163182441 


h(n) 


1 


0.810810810810811 


0.562063147746445 


0.371379608747416 


0.238302798822389 



The ratios a 2 (n) an d ore both strictly decreasing sequences with the limits 

lim a 2 (n) = lim (3 2 {n) = . (3.6) 

n— ►oo n— >oo 

Proof It is clear that a 2 (n) and /3 2 (n) cannot be negative. By Eqs. (13.11) - (13.31) . we 
have 

/ 2 (n + l) 



fl{n)g 2 {n) 
g 2 (n + l) a 2 {n 



6 + a 2 (n), (3.7) 
+ 7 + a 2 (n), (3.8) 



/ 2 (n)gf(n) p 2 {n) 

ho(n + 1) a 2 (n) / \ , 

\, \ = + 14 + 3ae,(n ■ 3 - 9 



Therefore, 



/ , ^ _ a 2 H^2H[6 + a 2 (ra)] _ g « 2 (n)[« 2 (n) + ff 2 (n)] 

a 2 (n) + 7f3 2 (n) + a 2 (n)/3 2 (ri) a 2 (n) + 7/3 2 (n) + a 2 (n)/3 2 (n) 

which shows that a 2 (n) is strictly decreasing. Similarly, we have 

a 2 (n)[a 2 (n) + 7/3 2 (n) + a 2 (n)/3 2 (n)] a 2 (n) a 2 (n)[3a 2 (n) - 7/? 2 (n)/3] 



12a 2 (n) + 14/3 2 (n) + 3a 2 (n)/3 2 (n) 3 12a 2 (ra) + 14/3 2 (n) + 3a 2 (n)/3 2 (n) 

(3.11) 
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With the initial values given in Table UTt 3/3 2 (n) < a2(n) is proved by induction. By Eqs. 
([3~TU]) and (jSXQ) , 

A(3 2 (n + 1) - a 2 (n + 1) 
_ a 2 {n) a 2 (n)[a 2 (n) + /3 2 (n)] a 2 (n)[12a 2 (n) - 28f3 2 (n)/3] 



3 a 2 (n) + 7(3 2 (n) +a 2 (n)(3 2 (n) 12a 2 (n) + U/3 2 (n) + 3a 2 (n)/3 2 (n) 

a 2 (n)X(n) 

' 3[a 2 (n) + 7(3 2 {n) + a 2 (n)f3 2 (n)]{l2a 2 (n) + Uf3 2 (n) + 3a 2 (n) (3 2 (n)] ' 

where 



(3.12) 



X(n) = l2al(n) + A%(3 2 {n)[l(3 2 {n)-a 2 {n)] + l2a 2 {n^ > 

(3.13) 

such that a 2 (n) < A/3 2 (n) is proved again by induction. Eq. (13.111) can be rewritten as 

R( ^\ = Rl \ Q2(n)[4^(n) ~ a 2 (n)][l + 3&(n)/4] + a 2 (n)(3 2 (n)[l - q 2 (n)/4] + Uf3 2 2 (n) 
) fj 2[ n) i2 ^ n) + u ^ + 3a2(n)/?2(n) 

(3.14) 

which shows that (3 2 {n) is strictly decreasing since a 2 (n) is less than four, i.e. its initial 
value, for all n > 1. Eq. ( 13.101) can be rewritten as 

<*(» + 1) = [l - ^MLzfM^ 1 , (3.15) 

7 L a 2 [n) + 7p 2 {n) + a 2 (n)p 2 (n)J 

which is always less than 6a 2 (n)/7 since fi 2 {n) is less than one, i.e. its initial value, for all 

n > 1 such that lim^oo a 2 (n) is zero. Finally, since (3 2 {n) < a 2 (n)/3, lim^oo (3 2 {n) is zero, 

and the proof is completed. □ 

We notice that the convergences of a 2 (n) and (3 2 {n) to zero as n increases are not rapid. 
The inequality 3f3 2 (n) < a 2 (n) can be improved a bit, and we state it as the following lemma. 



Lemma III. 3 For any n > no, 



3/3 2 (n) + ^<a 2 (n) , (3.16) 



where n = inf{n : a 2 {n) < 3/4, (3 2 {n) < 1/4} = 4. 

Proof By Eqs. (13TT0D and (13TTTD . we have 

„ / s otlin + 1) 
a 2 (n + 1) - 3f3 2 (n + 1) - 21 ^ ' 
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a 2 (n)[9a 2 (n) - 7/3 2 (n)] 



a2(n)[a 2 (n) + /3 2 (n)] 



12a 2 (n) + 14^ 2 (n) + 3a 2 (n)f3 2 (n) a 2 (n) + 7(3 2 (n) + a 2 (n)(3 2 (n) 

a 2 {n)(3 2 {n)[Q + a 2 {n)] 2 
27[a 2 (n) + 7p 2 (n) + a 2 (n)p 2 (n)} 2 

a 2 (n)Y(n) 

H»j + 7&(n) + a 2 (n)/? 2 (n)] 2 [12a 2 (n) + 14/? 2 (n) + 3a 2 (n)f3 2 (n)} 



where 



(3.17) 



K(n) = [a 2 (n) + 7/3 2 (n) + a 2 (n)(3 2 (n)}{ 3[a 2 {n) - 3/3 2 (n)][4^ 2 (n) - a 2 {n)} 



+3/3 2 (n) a 2 (n) - 3/3 2 (n) 
a 2 (n)p 2 (n) 



[3 + 2a 2 (n)] 



27 

'' [9a 2 2 (n) + 1008/3 2 2 (n) - 153a 2 (n) (3 2 (n) + Qa 3 2 (n) - 93a 2 2 (n) (3 2 (n) 



27 



-60a 2 (n)(3 2 (n) - Qa^n) (3 2 {n) - 50a 2 2 (n)f3' 2 (n) - 3ai(n)^(n)] 



> 3a 2 (n)[a 2 (n) - 3ft (n)] [3 + 2a 2 (n)] + 252ft (n)[4/3 2 (n) - a 2 (n)] 



27 



+a 2 (n)ft(n)[126 — 75a 2 (n) — 6al(n) 



60ft (n) - 50a 2 (n)ft(n) - 3c^(n)ft(n)]} . 

(3.18) 



Because Y(n) is positive whenever a 2 (n) < 3/4 and ft(n) < 1/4, which is true for all 
n > n = 4 by the previous lemma and Table [Til the inequality is established. □ 



We notice that although Eq. (I3.16P is by no means optimum, it is enough for the following 
lemma. 



Lemma III. 4 The sequence of the ratio {a 2 (n) / (3 2 {n)}^L decreases monotonically with the 
limit 

lim a 2 (n)/ft(n) = 3 . (3.19) 

Proof The initial value of the ratio is a 2 (0)/ft(0) = 4. It is clear from Eq. ( 13. lip that 
in the large n limit, the ratio a 2 (n)/f3 2 {n) is equal to three. By Eqs. (13.101) and ( 13. lip , we 



9 



have 



where 



a 2 {n) a 2 (n + 1) _ Z(n) 

~ fain + 1) ~ A(n)hW + + ' 



Z(n) = a 2 (n)[a 2 (n) + 7/3 2 (n) + a 2 (ri)/3 2 (n)] 5 



-/3 2 2 (n)[6 + a 2 (n)][12a 2 (n) + 14/? 2 (n) + 3a 2 (n)&(n)] 



(3.20) 



32 

[a 2 (n) - 3/3 2 (n)] 2a|(n)/3 2 (n) + —a 2 (n)(3%{n) + a 2 2 (n)(3%(n) 



+ 



a 2 (n) - 3/3 2 (n) 



a|(n) 
27 J 



[a|(n) + 17a 2 (n)/? 2 (n) + 2%f3 2 2 (n)} 



a 2 2 {n) 



o _ L a 2 2 (n) + I7a 2 (n)f3 2 (n) + 28^(n)] - -a 2 (n)^(n) . (3.21) 

With a 2 (n) > 3[3 2 (n), Z(n) is positive such that the sequence of the ratio decreases mono- 
tonically. □ 

The general expressions for / 2 (n) and g 2 (n) can be written as follows. 

Lemma III. 5 For a non-negative integer m and any positive integer n > m, 
/2H = / 2 (m) § ^(m) 2 jj |6 + a 2 (r 



n — 1) 



i=l 



n—m 

n[' 



x I I |7 + a 2 (n-j) + 
3=2 



a 2 {n- j) 



n—m 

3'" — 1 3'" — ""-)-i "■ — r r 

92(n) = f 2 (m) a ^(m) ~ j! 6 + a 2 ( 



(3.22) 



i=2 



n — z 



X 



n—m 

JJ [7 + a 2 (n-j) + 



Qi2{n- 3) 



(3.23) 



i/ere w/ien n — m = 1, the products with lower limit two are defined to be one. 



Proof It is clear from Eqs. (I3.7P and (13.81) that f 2 (m + 1) = f 2 (m)g 2 (m)[6 + a 2 (m)] 
and g 2 {m + 1) = / 2 (rn)sf (m)[7 + a 2 (m) + a 2 (m) / (3 2 (m)}. Consider Eqs. and fl3~23j) 

hold for a certain positive integer n = k, then 

f 2 (k + l) = f 2 {k)g 2 {k)[Q + a 2 {k)} 
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k — 771 



i=i 



k—m 



II 7 + a 2 (A:-i) + 



J=2 

x/ 2 (m) 

fc— m 

x JJ [7 + a 2 (A;-j) + 



Q2(fc - j) 



k—m 



~# 2 (m) 2 || |6 + Q! 2 (A; 

i=2 



«2(fc - j) 



[6 + a 2 (fc)] 



/ 2 (m) s ^ 2 (m) a [6 + a 2 (fc)][6 + a 2 (A; - l)] 1 



k—m 



i=2 
k—m 



3' + l 



J J 6 + a 2 (A;-z) 2 7 + a 2 (fc-l) 



a 2 (fc - 1) 



3=2 



11 W + a 2 (k-j) + 



®2(k - j) 

k—m+1 



3^-1 



3 fe — m + l^ 1 gfc — m + l —1 _ — _ r 

/ 2 (to) 2 ^ 2 ( m ) 3 || |6 + a 2 (fc + l 



i=i 



A:— m+1 



II [7 + « 2 (A; + l-j) 



i=2 



a 2 (A; + 1 - j) 
/5 2 (fc + 1 - j) 



(3.24) 



and Eq. (13.221) is proved by induction. Eq. (13.231) can be established by the same procedure. 
□ 

From the above lemmas, we have the following bounds for the asymptotic growth con- 
stant. 



Lemma III. 6 The asymptotic growth constant for the number of connected spanning sub- 
graphs on SG 2 (n) is bounded: 



\n[f 2 (m)g 2 (m)} + |ln60 
where m is a positive integer. 



< z S g 2 < 



Hf 2 (m)g 2 (m)] + \ ln[6 + a 2 (m)} h + a 2 {m) + 



(3.25) 
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Proof By Lemma 1111.5} we have 



<yn—m _|_ j gra— m j 

In /b^) = In fiijn) H m <7 2 (to) + A(n, to) 



(3.26) 



where 



A(n, m) = 



+ 1 



ln[6+a 2 (^-«)] + 5Z 9 ^ n 7+a 2 (™-j)4 



8=1 



J=2 



. (3.27) 



We have shown that as to increases, a 2 (fra) decreases to zero in Lemma IIII.2I and 
0:2(771) f fa (to) decreases to three in Lemma [III.4I such that 



A(n,m) < 

i=l 



"3' 1 • 1 



1 /3 n - rn - 1 



ri-m 3i _! _ l 

ln[6 + a 2 (m)] + ^ In 

J'=2 

+ n — m ) ln[6 + q^to)] 



7 + a 2 (m) + 



«2(mj 

AM 



1 /3 n ~ m - 3 
+ 2 ^ g n + m + 1 ) In 



7 + a 2 (m) 



a2\m) 

AM J 



(3.28) 



and 



A(n, rn) > 



ln6+ — - lnl ° 



1 / 3" - I 
2 



J'=2 

+ n- rnj In 6 + -( ^ n + TO+l)lnlO. (3.29) 



With the definition for zsg 2 given in Eq. (12.11) and the number of vertices of SG2(n) is 
3(3 n + l)/2 by Eq. (Q, the proof is completed. □ 



As to increases, the difference between the upper and lower bounds in Eq. (13.251) becomes 
small but the convergence is not rapid. We calculate the number of connected spanning 
subgraphs /2(to) up to to = 15, and we have the following proposition. 



Proposition III.l The asymptotic growth constant for the number of connected spanning 
subgraphs on the two-dimensional Sierpinski gasket SG2(n) in the large n limit is z,sg 2 = 
1.276495930.... 
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Without going into details, we state here without proof that the bounds can be improved. 
For a non- negative integer m and any positive integer n > m, the tighter bounds for 0:2(71) 
are 

d n - rn {m) < a 2 {n) < c n ' m (m) , (3.30) 

where 



c(m) 



6 + a 2 (m) 
10 + a 2 (m) 



d{m) 



6/3 2 (m) 



a 2 (m) + 7/3 2 (m) 



(3.31) 



It can be shown that 



{\n[f 2 (m)g 2 (m)] + -\n 



a 2 (m)i a 2 (m) rl 
42 + 6 *) { +t — v , ; X1 - + 



(5 2 (m) J [3 - c(m)\ lQ 7f3 2 (m) + a 2 (m) 

(3.32) 



and 



Z S G 2 > g^jH^Cm)^: 



In 60 Aa 2 (m) 
2\m)} + — — + 



17o|(m) 



2 15[3-d(m)] 900[3-rf 2 (m)] 
so that the asymptotic growth constant is zsg 2 = 1-27649593067.... 



} , (3.33) 



IV. THE NUMBER OF CONNECTED SPANNING SUBGRAPHS ON SG 2fi {n) 
WITH b = 3, 4 

The method given in the previous section can be applied to the number of connected 
spanning subgraphs on SGd,b( n ) with larger values of d and b. The number of configurations 
to be considered increases as d and b increase, and the recursion relations must be derived 
individually for each d and b. In this section, we consider the generalized two-dimensional 
Sierpinski gasket SG 2 ^{n) with the number of layers b equal to three and four. For SG 2)3 (n), 
the numbers of edges and vertices are given by 

e(SG 2 , 3 (n)) =3x6", (4.1) 
7 x 6" -I- 8 

v(SG 2 , 3 (n)) = , (4.2) 

where the three outmost vertices have degree two. There are (6 n — l)/5 vertices of SG 2j3 (n) 
with degree six and 6(6™ — l)/5 vertices with degree four. The initial values for the number 
of connected spanning subgraphs are the same as for SG 2 . f 2 ^(0) = 4, #2,3(0) = 1 and 
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h 2 , 3 (0) = 1. By the method illustrated in the previous section, we obtain the following 
recursion relations for any non-negative integer n. 

f 2)3 (n + 1) = / 2 6 » + lbfl 3 (n)g 2>3 (n) + 3/ 2 5 >)/i 2)3 (n) + 78/ 4 3 (n)4 3 (n) 

+18fl 3 (n)g 2 - 3 (n)h 2t3 (n) + U2f 23 {n)gl 3 {n) , (4.3) 

#2,30 + 1) = / 2 5 3 W#2,30) + f 2 , 3 (n)h 2i3 (n) + 16/ 2 4 3 (n)^ 23 (n) + 18f 23 (n)g 2t3 (n)h 2 ^(n) 
+89fl 3 (n)gl 3 (n) + 2/ 2 4 3 (n)^ 3 (n) + 77 fl 3 {n)gl 3 {n)h 2 , 3 {n) 
+171fl 3 (n)gi 3 (n) , (4.4) 

h 2 , 3 (n + 1) = 3/ 2 4 3 (n)^ 23 (n) + 6f 2t3 (n)g 2i3 (n)h 2 , 3 (n) + 51/ 2 3 3 (ra)^ 3 (ra) + 3/ 2 4 3 (n)/i 23 (n) 
+I29fl 3 (n)gl 3 (n)h 2 , 3 (n) + 279/ 2 2 3 (n)( 7 4 3 (n) + Q0fl 3 (n)g 2 , 3 (n)hl 3 (n) 

+564fl 3 (n)gl 3 (n)h 2>3 (n) + 468/ 2 , 3 (n)^ 5 3 (n) . (4.5) 

The figures for these configurations are too many to be shown here. Some values of f 2 , 3 {n), 
92,3{n), h 2t3 (n) are listed in Table HTT1 These numbers grow exponentially, and do not have 
simple integer factorizations. 



TABLE III: The first few values of f 2y3 (n), g 2l3 (n), h 2f3 (n). 



n 





1 


2 


/2,3H 


4 


56,192 


1,292,237,078,102,059,106,775,347,494,912 


52,3 (n) 


1 


24,624 


1,015,755,670,321,368,497,188,308,516,864 


h 2 , 3 (n) 


1 


33,792 


2,465,934,182,960,517,405,173,530,755,072 



The sequences of the ratios {a2,3(^)}.^Li and {P 2>3 {n)}'^ =1 defined in Eq. ( 13.41) again de- 
crease monotonically with lim^oo a 2)3 (n) = /3 2j3 (n) = 0. The ratio o; 2j3 (n)//3 2i3 (n) decreases 
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TABLE IV: The first few values of 02,3(71), /3 2) 3(n). The last digits given are rounded off. 



n 





1 


2 


3 


4 


"2,3 (jl) 


4 


2.28200129954516 


1.27219282732945 


0.660858801678112 


0.326587785819904 




1 


0.728693181818182 


0.411915158701392 


0.215917449918629 


0.107573237878269 



from four to three, the same as the results for SG^n). The values of 0:2,3 (ft-), ^2,3 (n) for 
small n are listed in Table IIVI 

By the same method as in Lemma IIII.51 we have the general expression for the number 
of connected spanning subgraphs. 

n—m 



/ 2 » = Wf^^m^n^n-j) 



(4x6 i_2 +l) 
1=1 



n^n-i)^" 1 - 1 ), (4.6) 

J=2 



where 



«9^(m) aL(m) „ „ 

P 2 , 3 (m) = 142 + + 78a 2 , 3 (m) + + 15a* 3 (m) + a| 3 (m) , (4.7) 

P2, 3 (^J P2,3(^J 

«9^fm) aL(m) aiJm) „ aL(m) 

W = 171+77f^+2-|^+89a 2 , 3 (m)+18-^+16a , 3 m +7^^+^ ,3 H 

(4.8) 

By the same argument given in Lemma IIII.61 we have the upper and lower bounds of the 
asymptotic growth constant for the number of connected spanning subgraphs on SG2${n): 

1 ■ [2 In / 2)3 (m) + 3 In g 2 , 3 (m)] + - 1 Rm [2 In 196 + 3 In 420] < z SG ^ 



7x gmL , n 35 x gr 

< -^[21n/ 2 ,3(m)+31n^3(m)] + ^^[21nP 2 ,3(m)+31ng 2 , 3 (m)] , (4.9) 
with m a positive integer. We have the following proposition. 

Proposition IV. 1 The asymptotic growth constant for the number of connected spanning 
subgraphs on the two-dimensional Sierpinski gasket SG2${ri) in the large n limit is zsg 2 3 = 
1.3972789680.... 
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We notice that the convergence of the upper and lower bounds remains slow. 



For S , G 2i4 (n), the numbers of edges and vertices are given by 

e(SG 2A (n)) = 3 x 10 n , (4.10) 

far, ,n 4xl0 " + 5 

v(SG 2A (n)) = , (4.11) 

where again the three outmost vertices have degree two. There are (10 n — l)/3 vertices 
of SG2,4(n) with degree six, and (10 n — 1) vertices with degree four. The initial values 
for the number of connected spanning subgraphs are the same as for SG2' /m(0) = 4, 
#2,4(0) = 1 and ^2,4(0) = 1. We wrote a computer program to obtain the recursion relations 
for SG2 A (n). They are lengthy and given in the appendix. Some values of f2,4(n), #2,4 ( n ), 
^2,4 (n) are listed in Table PVT. These numbers grow exponentially, and do not have simple 
integer factorizations. 

TABLE V: The first few values of f 2 iin), 52,4 ( n ), ^2,4 ( n )- 



n 


1 


2 




164,119,040 


27,140,375,625,882,898,681,725,275,604,427,985,839,201,951,967,246,962,831,668,354,270,630,763,784,348,631,040,000 


924 (n) 


77,622,016 


25,675,411,803,142,714,297,950,351,525,972,498,833,548,895,007,181,465,816,231,861,389,426,797,930,493,247,488,000 




112,848,672 


74,273,341,808,825,211,957,637,724,253,224,196,638,029,720,486,058,269,503,940,976,372,670,504,798,196,334,592,000 



The sequences of the ratios {o2,4(^)}£ ( Li and {P2,4( n )}n > =i defined in Eq. (13.41) again de- 
crease monotonically with lim^oo a 2 ^{n) = /5 2i 4 (n) = 0. The ratio 0:2,4 /Am M decreases 
from four to three. The values of 0:2,4 (w), A,4( n ) for small n are listed in Table [VTl 

By the same method as in Lemma IIII.51 we have the general expression for the number 
of connected spanning subgraphs. 

n—m n—m 

/ 2 , 4 (n) = / M (m)K^%( m) l(^) J] P 2 ,4(n-# (1 ° l ~ 1+2) II QzAn-j)^- 1 ^ , 

i=l j=2 

(4.12) 
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TABLE VI: The first few values of 02,4 (n), fo^in). The last digits given are rounded off. 



n 





1 


2 


3 


4 


a 2 , 4 (n) 


4 


2.11433622131123 


1.05705707211134 


0.475214294459902 


0.199993476915309 


/?2,4(n) 


1 


0.687841643364664 


0.345688118749653 


0.156412715166630 


0.662275943767262 



where 

P 24 ( m ) = H354 



al Am) / a\Am) 

■> , 1 2 #TT + 13626a 2j4 (m + 4140- 2 ^ 

*2,4 ( m ) $,4 ( m ) /VM 

«2 4( m ) 2 

:— — r + 6936a 24 
/3 2 4 (m) 



5856 + 516 ; 

/? 2 , 4 (m) P2,4l m , 

+174 ' ' 



^ 2 3 4 (m) Ti ™ 2 ^"^™/3 2 ,4(m) 

; 2 4 m + 1140 2 ' 4 ; ; + 15 J' 4 ; + 1928a 3 , 

ft,iH /3 24 (™) 

a« 4 (m) , (4.13) 



a^AJA 

a 24 (m) 
+144 2 ' 4 ( + 309a 4 
^2,4 (m) 



2 ! 4 (m) 

4, 4 (™) + 7^44 + 27« 2 5 4 (m) + a\ 



Q 2)4 (m) = 13732 + 14480 2 ' 4 ) ( + 2786 2 ' 4 ) + 82 J' 4 ) [ + 16250a 2 , 4 

/M m ) /? 2 , 4 ( m ) $2A m ) 

I 10609 " 2 ' 4(m) + 1095 + 12 a U m ' 

/3 2)4 (m) /5l i4 (^) /?|, 4 ( m ) 

a 4 



m 



9 / \ «2 4( m ) 

" L(m) + 3130 /SR 

aLM o , , aiAm) a\Am) 

+ 142 J' 4 + 2148a 3 4 m + 462 2 ' 4 + 6 -| 4 + 332a 4 4 m 
/? 2 4 (m) 2 ' 4V ^ 2i4 (m) /? 2 2 4 (m) 

-v 6 



K 24l m J 9 

r^— : — - + 8015a 2/| 
? 3 4 (m) 

m) 



«o 4 (m) , „ x an Am) 

+34-^-4 + 28a j4 m + + a , 4 m . 

/3 2 , 4 (m) P2,A m ) 

By the same argument given in Lemma IIIL6[ we have the 
asymptotic growth constant for the number of 

^-i^pn/a^m) + 2 In £ 2>4 (m)] + 36 x \ Q J kl 33620 + 2 In 84460] < z 5 g 2 , 4 



(4.14) 

21 we have the upper and lower bounds of the 
of connected spanning subgraphs on S , G 2)4 (n): 

f 2, 4 (^) + 2 In 5( 2 ,4 



jgraphs on SG2,±(n): 

[In 33620 + : 

36 x 10 m L 



with m 



a 



i/ 2 , 4 (m)+21ng 2|4 (m)] + 36 ^ — [In P 2 , 4 (m) + 2 In Q 2 , 4 
positive integer. We have the following proposition. 



[In f 2 4 1 

4xlO mL ^ ' 



m)] , (4.15) 



Proposition IV. 2 The asymptotic growth constant for the number of connected spanning 
bgraphs on the two-dimensional Sierpinski gasket SG2^(n) in the large n limit is zsg 2 a = 



su 

1.4849112 
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Here the convergence of the upper and lower bounds is again slow. 

V. THE NUMBER OF CONNECTED SPANNING SUBGRAPHS ON SG d {n) 

WITH d = 3, 4 

In this section, we derive the asymptotic growth constant of connected spanning sub- 
graphs on SGd{n) with d = 3,4. For the three-dimensional Sierpinski gasket SG 3 (n), we 
use the following definitions. 

Definition V.l Consider the three-dimensional Sierpinski gasket SG 3 (n) at stage n. (i) 
Define fz{n) = NcssG{SG 3 {n)) as the number of connected spanning subgraphs, (ii) De- 
fine g 3 {n) as the number of spanning subgraphs with two components such that one certain 
outmost vertices belongs to one component and the other three outmost vertices belong to 
another component. (Hi) Define h 3 {n) as the number of spanning subgraphs with two com- 
ponents such that two certain outmost vertices belong to one component and the other two 
outmost vertices belong to another component, (iv) Define r 3 (n) as the number of span- 
ning subgraphs with three components such that two certain outmost vertices belong to one 
component and the other two outmost vertices separately belong to other components, (v) 
Define s 3 (n) as the number of spanning subgraphs with four components such that each of 
the outmost vertices belongs to a different component. 

The quantities /^(n), ffetVi), h 3 (n), r 3 (n) and s 3 (n) are illustrated in Fig. [7J where only 
the outmost vertices are shown. There are four different classes of connected subgraphs 
enumerated by g 3 (n), three classes enumerated by h^n), and six classes enumerated by 
r 3 (n). The initial values at stage zero are /3(C)) = 38, #3(0) = 4, h 3 (0) = 1, r 3 (0) = 1 and 
s 3 (0) = 1. 

The recursion relations are lengthy and given in the appendix. Some values of fs(n), 
<7 3 (n), h 3 (n), r 3 (n), s 3 (n) are listed in Table NTH These numbers grow exponentially, and 
do not have simple integer factorizations. 
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h(n) g 3 {n) h 3 (n) r 3 (n) s 3 (n) 



FIG. 7: Illustration for the spanning subgraphs fz{n), gz{n), h 3 (n), r 3 (n) and s 3 (n). The two outmost 
vertices at the ends of a solid line belong to one component, while the two outmost vertices at the ends of 
a dot line belong to separated components. 

TABLE VII: The first few values of f 3 (n), g 3 {n), h 3 {n), r 3 (n), s 3 (n). 



n 





1 


2 


Hn) 


38 


8,554,560 


25,988,410,915,610,195,960,527,441,920 




4 


1,271,416 


4,544,490,996,892,396,578,747,598,336 


h 3 (n) 


1 


39,502 


73,629,059,909,730,939,289,401,696 


r 3 (n) 


1 


254,462 


917,115,147,969,863,922,701,973,216 


S3(n) 


1 


153,824 


637,427,406,318,067,141,227,862,784 



We find it is difficult to derive the bounds of the asymptotic growth constant for the 
number of connected spanning subgraphs on SG 3 (n). We calculate f 3 (m) up to m = 10, and 
fit the numerical value of the asymptotic growth constant to have the following proposition. 

Proposition V.l The asymptotic growth constant for the number of connected spanning 
subgraphs on the three-dimensional Sierpinski gasket SG 3 (n) in the large n limit is zsg 3 = 
2.06371. 

For the four-dimensional Sierpinski gasket SG^n), we use the following definitions. 

Definition V.2 Consider the four- dimensional Sierpinski gasket SG^n) at stage n. (i) 
Define f\{n) = NcssciSG^n)) as the number of connected spanning subgraphs, (ii) Define 
^ 4 (n) as the number of spanning subgraphs with two components such that two certain out- 
most vertices belong to one component and the other three outmost vertices belong to another 
component. (Hi) Define g' A {n) as the number of spanning subgraphs with two components such 
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that one certain outmost vertices belong to one component and the other four outmost ver- 
tices belong to another component, (iv) Define h^in) as the number of spanning subgraphs 
with three components such that one certain outmost vertices belong to one component, two 
certain other outmost vertices belong to another component and the remaining two outmost 
vertices belong to a third component, (v) Define h' A (n) as the number of spanning subgraphs 
with three components such that three certain outmost vertices belong to one component and 
the other two outmost vertices separately belong to other components, (vi) Define r±(n) as 
the number of spanning subgraphs with four components such that two certain outmost ver- 
tices belong to one component and the other three outmost vertices separately belong to other 
components, (vii) Define s A (n) as the number of spanning subgraphs with five components 
such that each of the outmost vertices belongs to a different component. 

The quantities fi(n), g^n), g'±(n), h^in), h' 4 (n), r 4 (n) and S4(n) are illustrated in Fig. 
[HI where only the outmost vertices are shown. There are ten different classes of connected 
subgraphs enumerated by ^(n), five classes enumerated by <? 4 (n), fifteen classes enumerated 
by h±{n), ten classes enumerated by h' A (n) and ten classes enumerated by r^in). The initial 
values at stage zero are / 4 (0) = 728, g A {0) = 4, ^(0) = 38, h 4 (0) = 1, h' A (0) = 4, r 4 (0) = 1 
and 54(0) = 1. 




h' 4 (n) r 4 (n) s 4 (n) 

FIG. 8: Illustration for the spanning subgraphs fi(n), gi{n), g'^in), h^in), h'^(n), r^n) and Si(n). The two 
outmost vertices at the ends of a solid line belong to one component, while the two outmost vertices at the 
ends of a dot line belong to separated components. 

We wrote a computer program to obtain the recursion relations for SG^n). They are 
too lengthy to be included here, and are available from the authors on request. Some values 
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of f±(n), g4(n), g'Jji), h^(n), h' 4 (n), r 4 (n), s 4 (n) are listed in Table IVTHl These numbers 
grow exponentially, and do not have simple integer factorizations. 



TABLE VIII: The first few values of / 4 (ra), g 4 (n), g 4 (n), /i 4 (n), /i 4 (n), r 4 (n), s 4 (n). 



n 


1 


2 


U(n) 


778,626,762,895,872 


1,024,406,418,765,003,907,906,096,145,114,250,200,136,082,865,744,856,739,402,552,777,712,697,606,144 


9i{n) 


88,489,486,528 


3,197,766,124,028,071,576,597,031,293,816,293,624,011,891,902,039,451,422,566,691,192,700,928 


9i( n ) 


52,683,007,497,792 


69,863,645,008,967,428,965,504,302,095,638,435,727,081,373,061,446,549,243,263,724,682,869,735,424 


h±(n) 


15,626,116,736 


482,200,982,250,980,661,780,613,757,386,400,225,524,299,614,671,798,975,981,818,835,107,840 


h' 4 (n) 


3,629,303,504,832 


4,765,691,494,696,738,414,627,738,223,389,422,884,987,040,021,355,690,863,179,839,765,575,892,992 


r 4 (n) 


258,767,297,696 


325,229,810,040,355,155,302,761,176,752,191,820,409,762,202,792,611,492,002,077,131,210,227,712 


54(77.) 


94,459,269,024 


110,974,534,976,153,854,286,043,758,382,592,092,762,465,813,295,695,669,459,295,951,908,765,696 



It is even more difficult to derive the bounds of the asymptotic growth constant for the 
number of connected spanning subgraphs on SG^(n). We calculate f^{m) up to m = 6, 
and we are satisfied with numerical fitting of the asymptotic growth constant to have the 
following proposition. 

Proposition V.2 The asymptotic growth constant for the number of connected spanning 
subgraphs on the four- dimensional Sierpinski gasket SG^n) in the large n limit is zsg 4 = 
2.7686. 



VI. DISCUSSION 



Compared with the asymptot ic growth constant for the number of spanning forests N$f 



on the Sierpinski gasket in Ref. 



Chang and Chen! ], we find that N C ssg is larger than N$f 



for all the considered cases. We conjecture that this inequality holds for all the generalized 
Sierpinski gasket. Define 

- Zo= lim '"^j G > . (6.!) 



V{G) 
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TABLE IX: Numerical values of z S q , zg, z SG d , and their ratios. The last digits given are rounded 
off. 



d 


D 


ZSG d 


ZG 


z SG d 


z SGd /z S G d 


ZSG d /zsG d 


2 


1.585 


1.048594857 


1.247337199 


1.276495931 


0.8214635326 


0.9771572077 


3 


2 


1.569396409 


1.666806281 


2.06371 


0.760475 


0.807675 


4 


2.322 


1.914853265 


1.981017076 


2.7686 


0.69163 


0.71553 



We list the first few values of ZsG d -> z SG d , an d their ratios in Table HXl 

As the spanning tree is a special case of connected spanning subgraphs where there is no 
cycles allowed, the number of spanning trees Nst{G) is always less than Ncssg(G). Define 

z_g = hm — — , (6.2) 

v(G)-*oa v(G) 

then z_ G < zg- We have obtained such asymptotic growth constants for the number of 
spanning trees on the Sierpinski gasket SGd for general d and SG2M with b = 3,4 in Ref. 



14j . They serve as the lower bounds for our current consideration for the connected spanning 



subgraphs. We list the first few values of z SGd , zsc d , and their ratios in Table HXl Notice 
that lower bound z SG is closer to the e xact value z<jg d w hen d is small, in contrast to the 



results for the spanning forests given in 
is large. 



Chang and Chenl ] that z SGd is closer to zq when d 
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APPENDIX A: RECURSION RELATIONS FOR 5G 2 , 4 (n) 

We give the recursion relations for the generalized two-dimensional Sierpinski gasket 
SG2,i{n) here. Since the subscript is (d, b) = (2,4) for all the quantities throughout this 
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section, we will use the simplified notation f n+ i to denote j '2,4 + 1) and similar notations 
for other quantities. For any non-negative integer n, we have 



f n+1 = C + 27f n g n + lf n h n + ?>mf n g 2 n + UAf n g n h n + 1928/ 7 ^ + 15/X + HW^X 
+6936/„ 6 «? 4 + inf n g n h 2 n + + 13626/ n 5 ^ + 2f n h\ + 516/ n 6 (?X 

+5856/ n ^X + H354/ n 4 ^, (Al) 

£n + i = /„V + /X + 28/ n 8 ^ + Uf n g n h n + 332/ 7 ^ + 6 /X + m2 f n g 2 n h n + 2148/ n 6 ^ 
+142/ n V^ + 3130 f n g 3 n h n + 8015/ n 5 ^ + 12/X + + 10609 f^X 

+ 16250/^ + 82 f n9n hl + 2786fSg 3 n h 2 n + U480f*g 5 n h n + 13732/ n 3 # 7 , (A2) 

h n+1 = 3f n g 2 n + Qf n9n h n + 87 f n gl + 3/X + 189/ 7 <?X + 1068/ n 6 «? 4 + 117 f n g n h 2 n 

+2558f^h n + UUflgl + + 18mf n g 2 n h 2 n + 17763/^X + 26934/^ 

+444/M + 12756/^X + 61422/ 4 «^ n + 53826/^ 7 + 20/X + 2388/ n 5 <?X 
+30948/ 4 sX + 83234/„ 3 ^„ + 42210/ 2 ^ . (A3) 
APPENDIX B: RECURSION RELATIONS FOR SG 3 (n) 

We give the recursion relations for the three-dimensional Sierpinski gasket SG^n) here. 
Since the subscript is d — 3 for all the quantities throughout this section, we will use the 
simplified notation / n+1 to denote f 3 (n + 1) and similar notations for other quantities. For 
any non-negative integer n, we have 

f n+1 = f n + \2f n g n + \2f n h n + \2f n r n + 48/ n 2 <? 2 + 9Qf 2 g n h n + 48/X + 72/„Vr„ 

+72f 2 h n r n + 56/ n <? 3 + lQ8f n g 2 n h n + lQ8f n g n h 2 n + 5Qf n h 3 n , (Bl) 
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9n+i = f n g n + 3/ n 3 r„ + 9f 2 g 2 n + 12 f n9n h n + f 3 s n + 36 f 2 g n r n + 30/^ n r n + 28/^ 
+66fn9n h n + ^fnQnK + 6/^ n s n + 6f 2 h n s n + 24/ 2 r 2 + 108/ n # 2 r n 
+ 192/ n ( 7 „/ in r ri + 84/ n /£ r „ + 20^ + 72<?X + 96<?X + 56 ^ n/l 3 ? (B2) 

/i n+1 = 2/ 2 /i 2 + 4/ n 2 /i n r n + 12/ n ^> n + 12/^Zi 2 + 16/X + 2/ n 2 r^ + 12/ n ^r n 

+48/ n <? n /i n r n + 36/ n /i 2 r n + 2g 4 n + 16g 3 h n + 36«?X + 32g n h 3 n + 22h 4 , (B3) 
r n+ i = f n gl + VlK + 6f 2 g n r n + 6/ n 2 /i n r n + 6f n g 3 n + 22 f n g 2 n h n + Uf n g n h 2 n + I6f n h 3 n 
+2f n g n s n + 2f n h n s n + 12/ 2 r 2 + mf n g 2 n r n + I32f n g n h n r n + 66 f n h 2 n r n + I2g 4 n 
+52^/i„ + 78<?X + A8g n h 3 n + 22h 4 n + 6f 2 n r n s n + Uf n g 2 s n + 28 f n g n h n s n 
+Uf n h 2 n s n + 120/ n <?„r 2 + 120/^r 2 + 88g 3 r n + 264g 2 n h n r n + 26Ag n h 2 n r n 
+S8h 3 n r n , (B4) 

s n +i = ±f n 9n + 36/ n # 2 r n + 2Af n g n h n r n + 12g 4 + 2Ag 3 n h n + I2f n g 2 n s n + 2Af n g n h n s n 
+I2f n h 2 n s n + U4f n g n r 2 n + I20f n h n r 2 n + U4g 3 r n + 360g 2 n h n r n + 216g n h 2 n r n 
+U4f n g n r n s n + \AAf n h n r n s n + 56g n s n + 168g n h n s n + 168g n h n s n + 56h n s n 
+208 f n r 3 n + 720g 2 n r 2 n + 1440^/i n r 2 + 720/i 2 r 2 . (B5) 
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